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Abstract 



The production of D* and D mesons in inelastic scattering of 160 GeV/c muons off a 6 LiD tar- 
get has been investigated with the COMPASS spectrometer at CERN for 0.003 (GeV/c) 2 < Q 2 < 
10 (GeV/c) 2 and 3 x 10~ 5 <x Bj <0.1. The study is based on 8100 events where a D° or 5° is de- 
tected subsequently to a D* + or D* decay, and on 34000 events, where only a D° or D° is detected. 
Kinematic distributions of D*, D and /T|(1430) are given as a function of their energy E, transverse 
momentum pj, energy fraction z, and of the virtual photon variables v, Q 2 and XBj- Semi-inclusive 
differential D* production cross-sections are compared with theoretical predictions for D* produc- 
tion via photon-gluon fusion into open charm. The total observed production cross-section for D** 
mesons with laboratory energies between 22 and 86 GeV is (1.9 ±0.4) nb. Significant cross-section 
asymmetries are observed between D* + and D* production for v < 40 GeV and z > 0.6. 
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1 Introduction 

The production of D mesons in inelastic scattering of 160 GeV/c muons on nucleons fiN — > jJ.'DX is 
assumed to be dominated by a process where the exchanged virtual photon y* fuses with a gluon g 
into a charm anti-charm quark pair, y*g — > cc. The cross-section a^ g ^ cc of this photon-gluon fusion 
(PGF) process and its dependence on the relative polarization of photon and gluon can be calculated in 
perturbative QCD jTHS]]. Thus, using polarized muons and polarized nucleons, a measurement of the 
photon nucleon cross-section asymmetry Aa^ N ^ ccX / G' N ^ ccX allows the determination of the gluon 
polarization Ag/g in the nucleon. With this objective, open charm production has been studied in the 
COMPASS experiment at CERN for longitudinally polarized muons interacting with longitudinally po- 
larized deuterons. The incoming muon energy of 160 GeV was chosen, since the cross-section difference 
AgY* n ^ ccX for parallel and anti-parallel spins of photon and nucleon reaches a maximum for virtual 
photon energies around 80 GeV according to most models for the gluon helicity distribution function 
Ag(xsj, Q 2 ), and the polarization transfer from muon to virtual photon is large in the relevant photon 
energy range. 

Final states, where the decays D° — > K~n + or D* + — > D Q n + — > K~K + K + or the charge conjugate decays 
are detected, are chosen in order to achieve the best possible combination of mass resolution, signal- 
over-background ratio and signal statistics. Based on samples of events with these final states, extracted 
from data taken during the years 2002-2006, COMPASS has published results for (Ag/g) 0. 

The photon-gluon cross-section asymmetry an = Ac 7 * g ^ cc / c^* 8 ^ cc needed for extracting (Ag/g) is 
estimated making two assumptions: only PGF contributes as calculated in leading order QCD and charm 
and anti-charm quarks hadronize independently of the target polarization. The parton kinematics are 
estimated event-by-event on the basis of the observed 3-momentum of the D° meson and the momentum 
difference of the incoming and the scattered muon using a parametrisation based on the Monte Carlo 
event generator AROMA with default charm quark fragmentation [10]. 

However, production mechanisms other than PGF with standard charm quark fragmentation may con- 
tribute to the observed events with charmed mesons. The interaction of the virtual photon with intrinsic 
charm lfTTT4T3Tl is one possible competing mechanism. The associated production of A C D° lTT4ll or, more 
generally, asymmetric hadronization of c and c like in the Dual Parton model with a meson and a baryon 
string lfI31[T6l may play an important role in some regions of phase space. A study of the phase space dis- 
tributions and semi-inclusive differential cross-sections of the D and D* mesons within the acceptance 
of the COMPASS spectrometer may yield information about the contributions of different production 
mechanisms ifTTl . 

At HERA, i.e. at much larger center-of-mass energy, charm electro-production has been studied in detail 
by the HI and ZEUS Collaborations, see IPT&I - I201 and references therein. In addition to the PGF other 
production mechanisms also contribute in this case, like gluon-gluon fusion to cc from a resolved pho- 
ton. The hadronizations of c and c can more safely be assumed to be independent. COMPASS covers 
a complementary kinematic region with virtual photon energies in the range from threshold at about 30 
GeV up to 140 GeV in the laboratory frame. Prior to COMPASS, this energy range was covered only by 
the EMC experiment ETTl . which collected about 90 D° meson events produced by deep inelastic scat- 
tering of 240 or 280 GeV muons on hydrogen and deuterium targets for a study of the charm production 
mechanism. Only one charm photo-production experiment explored the region close to threshold ll22l . 
while two concentrated on high energy photons |[T6ll23l . 

The present article shows details of the phase space distributions of D meson production as a function 
of various kinematic variables: the energy v = — EL of the exchanged virtual photon y* (assuming 
single photon exchange) with four momentum q = p — p', the inelasticity y = v/E^ of the event, the 
negative invariant y* mass squared Q = —q = —(p — p') and the Bjorken scaling variable xsj = 
Q 2 /2P ■ q = Q 2 /2M p v. Here E^ and E'^ are the laboratory energies, p and p' the 4-momentum vectors 
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2 EXPERIMENTAL SETUP 



of the incoming and scattered muon respectively, P is the 4-momentum of the target nucleon and M p is 
the proton mass. 

In order to describe both D* and D° meson production, the following kinematic variables are used: the 
transverse momentum pj of the Kn pair (from the D° decay) with respect to the y* direction, the D° 
energy E in the laboratory system and the energy fraction z = E/v. 

The outline of the paper is as follows: after a brief overview of the experimental set-up (section [2]) and 
a detailed description of the data selection procedure (section [3]>, the methods of signal extraction are 
described in section [4] The kinematic distributions of events from signal and background regions are 
shown in section [5] They are based on the entire available data sample collected during the years 2002- 
2006, and are not corrected for acceptance. The purpose of this section is to compare the distributions of 
open charm to those of background events. Section[6]describes the acceptance correction and luminosity 
calculation needed to extract the total and differential semi-inclusive cross sections for charm production. 
They are performed for the 2004 data only. In section[7J the differential cross-sections as a function of the 
various kinematic variables and the total cross-section obtained by integration are shown and compared 
with available theoretical (AROMA) predictions for the production of D mesons by PGF. Significant 
differences between D* + and D* meson production are observed for the acceptance-corrected data 
from 2004. A statistically more precise comparison of D* + and D* production is based on the entire 
data sample (2002-2006). Particle-antiparticle asymmetries are determined under the assumption, which 
is verified for the 2004 data, that the D* + and D* acceptances are equal to a good approximation. 

2 Experimental setup 

The data were taken using the COMPASS spectrometer situated at the M2 beam line at the CERN Super 
Proton Synchrotron. A detailed description of the COMPASS spectrometer can be found in Ref. [24]. 

The momentum of the positive muon beam is about 160 GeV/c with a spread of 5%. The momentum of 
each incoming muon is measured with a precision of Ap/p < 1% in the beam momentum spectrometer 
located upstream of the experimental hall, and its direction is measured with a precision of 30 /irad with 
a detector telescope in front of the target. 

The (polarized) 6 LiD target consisted of two 60 cm long cells during the years 2002-2004 and of 3 cells 
with a total length of 120 cm in 2006. The polarization is reversed regularly such that the products of 
integrated luminosities times acceptance are equal for both polarizations. The sum of both corresponds 
to the unpolarized case. Hence unpolarized distributions, which are the subject of the present analysis, 
are obtained from the sum of all data. The target is housed in a superconducting solenoid magnet, which 
determines the angular acceptance of the spectrometer. The acceptance in the polar angle, measured at 
the upstream edge of the target, was 70 mrad in 2002-2004, while with the new target magnet in 2006 it 
was increased to 180 mrad. 

The 2-stage spectrometer is designed to reconstruct and identify the scattered muon and produced hadrons 
over a wide momentum range. It contains a large angle (LAS) and a small angle (SAS) part, each part 
equipped with a dipole magnet. Tracking detectors are located in front and behind each magnet, and 
electromagnetic and hadron calorimeters behind. The LAS covers polar laboratory angles from about 15 
mrad up to 70 mrad in 2002-2004 and, with the new target magnet, up to 180 mrad in 2006. The SAS 
covers the polar laboratory angles below 20 mrad. 

Particle tracking is done with a large variety of tracking detectors: several stations of silicon microstrip 
detectors, scintillating fiber detectors, high resolution micromesh gaseous chambers, gas electron-multi- 
plier chambers, drift chambers, large area straw drift chambers, multiwire proportional chambers and 
muon drift tubes. The scattered muons are identified downstream of additional hadron absorbers placed 
behind the hadron calorimeters. Charged hadrons are identified by a Ring Imaging Cherenkov detector 
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(RICH) in the LAS. 

The trigger system [25 ] uses hodoscope and calorimeter information to select inelastic muon interactions 
with minimum bias. The overall trigger and muon track reconstruction efficiency is in the range 60% to 
80% for most of the kinematic region covered by COMPASS. 

3 Data selection 

The total number of events with an incoming muon (140 GeV/c < < 180 GeV/c) and a scattered 
muon from a common vertex is 5.2 x 10 9 , which corresponds to an integrated luminosity of about 2.8 
fb _1 . This sample is used to search for D°(D°) and D** mesons. A fiducial volume cut makes sure that 
the extrapolated incoming muon trajectory traverses all target cells and that the primary vertex is located 
within the volume of one of the target cells. 

Since the COMPASS experiment uses a large solid target, the detection of a secondary decay vertex, 
which is a standard method in open charm detection, is excluded and the selection of D°(D°) and D** 
mesons relies on requirements on event kinematics and particle identification. The event selection is the 
same as used in the previous COMPASS open charm publication [9] except for stricter requirements on 
the selection of the incoming muon. 

Cuts used to select D° originating from the decay of a D* (referred to as D* or 'tagged' D° sample) 
slightly differ from those used to select directly produced D° mesons (referred to as 'untagged' D° 
sample). An event from the untagged D° sample contains at least one candidate for the 2-body decay 
D° — > K~7i + or its charge conjugate (c.c), while in the tagged D sample a slow pion from the decay 
chain D* + — > D°7r+ — > K~n + nf (or c.c) has to be present in addition. 

Particles are identified by using the RICH. All tracks with momentum measured in one or both spec- 
trometer stages and falling within the geometrical acceptance of the RICH are used to calculate the 
likelihoods (LKs) that the Cherenkov photons detected by the RICH are due to electron, muon, pion, 
kaon, proton, or background. The LK for a specific particle is calculated only if the particle velocity is 
above the threshold for the emission of Cherenkov photons in the radiator gas. This threshold depends on 
the refractive index that is extracted from the data on a run-by-run basis. For pions, kaons, and protons, 
this gives an average momentum threshold of 2.5, 8.9 and 16.9 GeV/c respectively. At large momenta 
pions and kaons cannot be efficiently separated, thus it is required that the momentum of the particle 
is below 50 GeV/c. In the tagged D° sample, due to the small mass difference between the D*(2010) 
and the D°(1865), only a limited energy is available for the pion produced in the D* — > D Q n s decay. In 
this case, the n s candidate must not have been identified as an electron by the RICH. Details on the LK 
requirements and the use of the RICH information can be found in Ref. [ 17 ]. 

For untagged D°, the following cuts are applied to the K~7i + and K + 7i~ combinations: p n >l GeV/c, 
z > 0.2, \M(Kn) -M(D°)\ < 700 and \cosd K \ < 0.65, where B K is the decay kaon angle in the D° 
center-of-mass system with respect to the D° direction of flight. 

For the tagged D°, the Knn s invariant mass is calculated only if the K% system has an invariant mass 
in the range \M(Kri) —M(D°)\ < 700 MeV. The distribution of AM = M{Knn s ) - M(Kn) -M(n) as 
a function of M(Kn) is shown in Fig. [I] Here a clear spot for the D* is visible at AM ~ 6 MeV in the 
region of the D° mass. The cut 3.2 MeV < AM < 8.9 MeV improves the D° signal with respect to the 
combinatorial background by more than an order of magnitude. The Kn system is also required to have 
z> 0.2 and | cos 0^1 < 0.9. 

These sets of cuts define the untagged and tagged D° samples, i.e. the D° and D* candidates. 

The cos 6k distribution is the only distribution where a safe theoretical prediction can be made. The 
uncorrected cos 6k distribution of K% events before any mass cuts, i.e. mostly background, shown in 
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100 



50 



1.4 



1.6 



2 2.2 2.4 

M(Kn) (GeV) 



Fig. 1: Scatter plot for D* candidates before applying the AM cut. Vertical axis: AM; horizontal axis: 
M(Kn). The accumulation of events around the D nominal mass of 1.864 GeV and AM = 6.1 MeV 
corresponds to the decay sequence D* — > n s D° — > % S (K%). 
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Fig. 2: Distribution of cos(fc in the Kn rest frame for (mostly) background combinations (scaled by 
0.001, solid line), for the D° signal region (scaled by 5, full circles) and for the acceptance corrected D° 
signal (open squares). The dashed lines correspond to the | cos 6k\ < 0.9 cut. The D signal is from the 
2004 tagged D° sample. 



Fig. [2] is strongly peaked towards cos 6k = — 1- For signal events, the cos0^ distribution should be flat 
after acceptance correction since the D° has spin 0. This expectation for the D° is confirmed in Fig. [2] 
where, for the tagged D° sample, the distribution for D° is shown before and after acceptance correction 
(the method of signal extraction and the correction for the acceptance will be described in sees. [4] and [6]). 

The so-called ambiguity cut applied in Ref. ||9l is not applied in the present analysis. This cut discards an 
event if two D or D meson candidates are found within the mass window of ±700 MeV and removes 
a significant number of good events. However, the probability to find two D° (or two D°) mesons in the 
signal peak is practically zero. Hence the present analysis, which extracts the number of signal events 
from fitting separately the D° and D° peaks, does not suffer from this ambiguity. 

In the mass window of ±700 MeV around the nominal D° mass, the tagged D sample consists of 160 x 
10 3 neutral Kn combinations. In order to avoid overlapping samples, at this stage, the D* candidates 
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xlO 3 xlO 3 




M(Kn) (GeV) M(Kn) (GeV) 

Fig. 3: Invariant M{Kn) mass spectra within a window of ±700 MeV around the nominal D° mass, (a) 
D° sample before and (b) after background subtraction, (c) D* sample before and (d) after background 
subtraction. Both neutral charge combinations are shown separately, together with their sum in (a) and 
(c). See text for the background subtraction by fits. 



are taken out of the untagged D sample. In the same mass window, the untagged D sample comprises 
17 x 10 6 neutral K% combinations. 

The invariant K% mass spectra are shown in Fig. [3^ for the untagged D° sample, for all neutral Kn 
combinations and also separately for the K + n~ and K~7i + combinations. These spectra exhibit the D° 
peak at 1865 MeV. The prominent peak to the left is due to the decay of the narrow ^|(1430). In Fig. [5]: 
invariant mass spectra are shown for the tagged D° sample. In this case, only some feed-through of the 
^(1430) resonance is seen and a pronounced, rather narrow peak about 250 MeV below the nominal 
D° mass. As shown by Monte Carlo simulations, this peak at 1620 MeV results from 3-body decays of 
the D° — > KnnP, where the n° escaped detection, with some contributions from D* decays with more 
than 3 particles in the final state. The signal-over-background (S/B) ratio is about 1:1 for the events of 
the tagged D° sample. For the untagged D° sample, S/B is only 1:10, but the number of signal events is 
four times higheiQ 

Mass spectra of all the Kn combinations are shown in Fig. [4^ separately, using only data from 2004. 
The spectra for the two neutral charge combinations show three narrow peaks corresponding to K*(S90), 
#2* (1430) and D°(1865). Also, other short lived kaonic (strange) resonances are present but they super- 
impose together with combinatorial background to a structureless distribution that can almost perfectly 
be described by a single exponential function, see Fig|4j). 



The S/B is calculated in a ±50 MeV mass window (±2cr) around the D° peak. 
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r 

0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 

M(Kn) (GeV) 

Fig. 4: Kk invariant mass distribution for the full mass range using data from 2004 only, (a) All four 
charge combinations (from top to bottom + -,-+, — , + +) are shown separately, with a linear vertical 
scale; (b) only neutral, non-exotic charge combinations (- +, + -) summed up, logarithmic vertical scale. 
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4 Method of signal extraction 



The invariant Kn mass spectra shown in Figs. [3^,c are fitted with a function given by the sum of the 
following elements: a Gaussian for the D° — > Kn signal, an exponential for the background, a shape 
determined by Monte Carlo simulations for the peak at 1620 MeV from 3-body decays of the D° and 
by relativistic Breit-Wigner intensities for the XJ(1430) an ^ ^3 (1780). The latter K resonance is barely 
visible in this spectrum but shows up clearly in certain kinematic regions, see below. 

The Kn spectra are remarkably well described fitting them with 12 parameters, as described above. 
Figs. [3]),d show the spectra after subtraction of the exponential background. From the fits one obtains 
(3610 ± 90) D*+ -> (K-7l + )nf and (4530 ± 100) D* -)• (K + n~)7V~ for the tagged sample as well as 
(15200 ± 800) D° -> K-%+ and (18400 ± 900) D° -> K + tc~ for the untagged D° sample. 

The dependence on kinematic variables of the production rate of D° and D*, together with those of the 
neighbouring A|(1430) resonance and the background, is extracted by fitting the mass spectra for each 
kinematic bin. Alternatively, the signal distributions of the D° and the D* are obtained by side-band 
subtraction. 

Using the first method, the fit yields in every bin of a given kinematic variable the number of D, K\ (1430) 
and (1780) together with the background. In Figs.[5^-f examples of the Kn invariant mass spectra for 
different intervals in z are shown before and after the subtraction of the fitted exponential background. 
These fits did not include the X|(1780). The fitting method allows monitoring of all details of the fit, as 
illustrated in the inserts of Fig. [5] The broad structure showing up for z > 0.75 in Fig. [5] is attributed to 
the 7^3 (1780) resonance. This resonance follows the same behaviour as the X|(1430) resonance, i.e. it 
is produced at larger values of z than the D° (see sec. [5]). The introduction of the K^ (1780) resonance in 
the final fit also removed a small but statistically significant and unexplained discrepancy between fit and 
data on the left side of the D° — > Kn peak in the z-integrated spectrum, where the fit before the inclusion 
of the X| (1780) was systematically below the data (see Fig. [5^ to be compared with Fig. J5J}, where the 
£3" (1780) has been included). 

The second method for signal extraction is the standard side-band subtraction. This method can only be 
applied to the D° and the D* signals, due to the limited mass range (±700 MeV around the nominal D° 
mass). Three Kn mass windows are chosen. The central one, which is 100 MeV wide and centered at 
the nominal D° mass, contains the D — > Kn signal plus background. The two side-bands contain only 
background events. They are 50 MeV wide and centered at ±100 MeV above or below the nominal D° 
mass. Thus three independent distributions are obtained as a function of each kinematic variable. The 
sum of the side-band distributions is subtracted from the central distribution, assuming that the side- 
band distributions correctly represent the distribution of background under the signal. This assumption 
is supported by the observed similar behaviour of the distributions in the two side-bands. 

Usually, the background below the signal is obtained by linear interpolation between the side-bands. 
Such a linear interpolation overestimates the background under the signal. Therefore it cannot be applied 
for the untagged D° sample, where S/B ~ 1/10. Instead, an estimate of the background under the signal 
is obtained from the fit. The total number of background events in the two side-bands is correspondingly 
rescaled. 

For the chosen width of the central window, about 5% of the signal is found outside. Hence the number 
of signal events obtained by side-band subtraction is expected to be lower by 5% than that obtained with 
the signal fitting method. 
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4 METHOD OF SIGNAL EXTRACTION 
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Fig. 5: a-f) Invariant Kn mass spectra in bins of the energy fraction z for the untagged D° sample. 
The vertical (red) lines indicate the nominal positions of ^(1430) and D°. The inserts demonstrate 
the signal behaviour after the removal of the fitted exponential background. The fit contains D° — > Kn 
at 1865 MeV, D° ->■ #7171° at 1620 MeV, the ZJ(1430) and an exponential background, g) shows the 
signal behaviour after removal of the fitted exponential background for the entire z range (no ^3 (1780) 
assumed). This figure has to be compared with Fig. |5Jd where the £3* (1780) was included in the fit. 
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5 Comparison of kinematic distributions 

In this section, event distributions are shown as a function of the relevant kinematic parameters, for both 
the tagged and untagged D° samples as well as for the X|(1430) and background. The data collected 
in 2002-2006 are used, and the distributions are not corrected for acceptance. However, the geometric 
acceptances for the various compared Kn systems are similar. 

The distributions of the X|(1430) signal are obtained from the untagged D° sample using the signal 
fitting method. The distributions of the Kn -background combinations are extracted from the two Kn side 
bands of the tagged D° sample, at invariant masses of 1765 ± 50 and 1965 ± 50 MeV. The kinematic 
distributions of D° and D* are obtained by applying both signal extraction methods described above, 
allowing to cross-check the stability of the result. While for the tagged D° sample perfect agreement is 
found between the two methods, for the untagged sample some disagreement beyond the statistical error 
is observed, for instance at low values of z or low Kn energy E. This is the result of strongly varying 
background shapes with additional broad resonances emerging below the XJ(1430). The corresponding 
data points for D° and K^ ( 1430) are omitted, since a more complex background description would be 
needed. 

In Fig. [6] the distributions of the D°, the ZJ(1430) and the background under the D° are compared, 
showing their different behaviour. The distributions as a function of the inclusive variables Q 2 and XBj 
are displayed in Figs. [6^,b. For the tagged sample, the average values of Q 2 and xbj extracted from 
these distributions are about 0.5 (GeV/c) 2 and 0.005, respectively. Some differences between signal 
and background events are observed at large values of Q 2 and xgj. As a function of v, the distributions 
for the various Kn systems are significantly different (see Fig.|6j:). The £|(1430) distribution peaks at 
lower values than that of D°, and the rise at low v that is caused by the increase of both acceptance and 
cross-section starts at lower v. The background peaks at a somewhat higher values, but has a similar rise 
with v as D° and D*. 

No clear differences are observed between the distributions as a function of the Kn energy E (see Fig.[6]l). 
Given the reason described above, the K\ ( 1430) and the D° data points at lower values of E are omitted 
from the untagged sample. The distribution for the D* signal as a function of p\ (Fig. |6^) shows an 
almost single-exponential decrease, while the distribution for the D° flattens above 3 (GeV/c) 2 . The 
difference between D° and D* may be related to the fact that for the D* only the p\ of the 2-body sub- 
system is shown. Both distributions are significantly different from those of background and A^(1430). 
From a fit of an exponential function up to p\ = 2(GeV/c) 2 , the following slopes are obtained in units 
of (GeV/c)- 2 : -0.84±0.03 forD*, -0.96 ±0.06 for D°, -1.94±0.01 for^ 2 *(1430) and -1.69 ±0.01 
for background. The distributions in z show significant differences, too. The background is concen- 
trated at smaller values of z than the D° signal. Moreover, the distribution of the X|(1430) is peaked at 
significantly higher values of z than that of the D°. 

In conclusion of the comparison: remarkable differences are observed between the distributions of the D 
meson signals, the iT|(1430), and the background as a function of the kinematic variables v, pr and z. 
This clearly points to different production mechanisms for D mesons and the £|(1430). The observed 
differences between D mesons and K^ ( 1430) agree qualitatively with the differences expected if the D 
mesons result from the fragmentation of a pair of charm quarks and the X| (1430) from the fragmentation 
of a quark knocked out in a leading order process. 

The interpretation of the kinematic distributions of the background is more complex, since this back- 
ground is dominated by combinatorial entries. No attempt is made to interpret it. However, one should 
mention that other background events of non-combinatorial origin (e.g. in the untagged sample the back- 
ground taken from side bands has also large contributions from resonances or from nK correlated pro- 
duction in the fragmentation) have been observed to behave very similar to the background shown in 
Fig-© 
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Fig. 6: Measured kinematic distributions of various (Kn) systems before acceptance correction as a 
function of (a) Q 2 , (b) x Bj , (c) v, (d) E, (e) p\, and (f) z. The symbols D° and K\ denote D° and 
( 1430) from the untagged sample. The symbols D* and bkg denote D* and background from the two 
side-band windows for the D* sample. The data are from the years 2002 to 2006. 



6 Acceptance and integrated luminosity 

Acceptances and integrated luminosity, which are needed to extract semi-inclusive total and differential 
cross-sections, are calculated only for the tagged D sample of the year 2004. Since this is the first 
detailed acceptance calculation for this particular final state at COMPASS, the present section also aims 
at illustrating the acceptances of the COMPASS spectrometer for the detection of the scattered muon and 
the D*^. For this reason, 2-dimensional acceptances will be shown as a function of selected variables. 

Acceptance calculations are done using a complete Monte Carlo simulation of the detector configuration, 
including the triggers and the track reconstruction code for the 2004 data. Events are generated using 
AROMA 2.2.4 ifTOl . which assumes photon-gluon fusion into cc to be the dominant underlying mecha- 
nism for D* production. Default fragmentation functions are used and parton showers are generated. The 
charm quark mass is set to 1.35 GeV. Produced D*s are forced to decay to D°7i + — > K~n + 7t + for D* + or 
to D°7t~ — > K + 7i~7i~ for D* . Trigger conditions and data selection criteria applied to the Monte Carlo 
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Fig. 7: Number of generated (AROMA) events as a function of (a) xsj and y and (b) p± and E. 
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Fig. 8: (a) 'Inclusive' acceptance A fl (xsj,y) and (b) overall acceptance Ad* (xBj,y) before applying thei? 
cut. 



events are the same as for real data. In total, 10 7 events were generated for both decays. The acceptances 
are calculated as a function of the reconstructed values of the kinematic variables, thus accounting for 
experimental resolution and bin-to-bin smearing. 

Figure [7] shows the number of generated events (a) as a function of XBj, y and (b) as a function of p±, E 
of the D meson. The transverse momentum p± is measured with respect to the direction of the incoming 
muon beam. In both pictures the generated events are mainly concentrated in the lower left corner. 

For illustration, the acceptance for D* production is shown at two stages, i.e. after requiring the recon- 
struction of the scattered muon and after the additional reconstruction of the three hadrons from the D* 
decay. The 'inclusive' acceptance A^(xBj,y) is shown in Fig. [8^, and the overall acceptance Ao*(xBj,y) 
in Fig. [8j). In the kinematic region relevant for charm production, the inclusive acceptance A^(xBj,y) 
is fairly homogeneous and ranges between 50% and 80%. The overall acceptance Ao*(xsy,j) is also 
homogeneous for y > 0.2 and ranges from 1% to 5%. The cut-off at y = 0.2 is due to the momentum 
selection for the RICH identification. 

The overall acceptance Ad* as a function of E and p± (i.e. transverse momentum with respect the in- 
coming muon) is shown in Fig. [9] The upper limit of about 100 mrad for the spectrometer acceptance 
in the year 2004 can be seen at low energy and large p±. For 20 GeV < E < 80 GeV the acceptance 
ranges between 5% and 13%. Outside this energy region the acceptance drops to zero due to the lack of 
particle identification and therefore 20 GeV < E < 80 GeV is required in the further analysis. The one- 
dimensional acceptances used below to determine the differential inclusive cross-sections are limited to 
this range of D° energies. 



The one-dimensional acceptance functions A D *+ and A D *- are shown in Fig. 10 as a function of v, E, z 
and pj. In addition, the ratio of the D* + acceptance over that of D* is shown in each case. Note that 
within the statistical uncertainty of the Monte Carlo simulation the acceptances are almost equal for D* + 
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Fig. 10: One-dimensional acceptances for D* 1 * 1 production as a function of (a) V, (b) (c) z and (d) 
/7y. Red boxes correspond to D* + , blue triangles to D* events. The black circles show the ratio of 
acceptances A D »+ /A D *-, the ordinates for the ratios are drawn on the right 4iand side of the figures. 



and D* , with the acceptance for D* + being slightly higher than that for D* . 

The integrated luminosity Jz? is determined by a comparison of the measured number of inclusive in- 
elastic muon scattering events with the best available measurement of the corresponding cross-section. 
The differential number of events is the product of integrated luminosity, inclusive muon acceptance and 
inclusive differential cross-section: 

*N(x BJ ,y) _ „ A d^ N ^(x Bj ,y) 
dxgj dy 1 ' dxgj dy 
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The inclusive inelastic muon-deuteron cross-section was measured by the NMC Collaboration for various 
muon energies between 90 and 280 GeV and published as a parameterization of the structure function 
Fi 11261 . Thus the cross-section has to be reconstructed based on this F2 parameterization. The measured 
cross-section is connected with the one -photon exchange cross-section via a radiative correction factor 

d 2 a^' x (x B j,y) _ 1 d 2 a ly (x Bj ,y) 



dx B jdy V( x Bj,y) dx Bj dy 

The one-photon exchange cross-section is connected with F2 by : 

d 2 a l7 (x Bj ,y) = An{ahcf ( Q 2 (l-2m 2 /Q 2 )(y 2 + Q 2 /E 2 ) \ 

dxBjdy x Bj yQ 2 I J 4E 2 + ' 2(1 +R(x Bj ,Q 2 )) J 2 ^'^ 



(2) 



(3) 



where m is the muon mass. The factor R(x B j, Q 2 ) is the cross-section ratio for longitudinal over transverse 
photons: 

R(xbj,Q 2 ) = — ■ (4) 

The radiative correction factor v\ (xBj,y) is calculated with codes based on |[27l . The ratios R(xgj, Q 2 ) are 
determined as in Ref. II281 . Given the light material composing the target (Li, D and He), nuclear effects 
have been neglected. 

The integrated luminosity is determined in bins of (x B j,y) as: 

a, = 1 d 2 N(xBj,y)/(dx Bj dy) 

~ A^x Bj ,y) ' d 2 a^' x (x Bj ,y)/(dxBj dy)' 1 } 

The integrated luminosity on the left-hand side of Eq. [5] has to be constant, while all terms on the right- 
hand side depend on x B ; and y. As a side product of extracting the integrated luminosity, this equation can 
be used to evaluate the uncertainty of the muon acceptance calculation for Q 2 values larger than about 
0.6 (GeV/c) 2 , where the NMC parameterization is valid. The values of Jz? obtained for different (xgj,y) 
bins vary indeed by up to 20% over the relevant (x B j,y) range, so that an overall systematic uncertainty 
of 20% is attributed to the product of integrated luminosity and inclusive muon acceptance. The average 
value of the integrated luminosity is calculated as a weighted mean of the luminosities determined in 
(xsj,y) bins, using the data at Q 2 > 0.6 (GeV/c) 2 . For a given bin the weight is the number of events in 
that bin. The result for the integrated luminosity of the 2004 data is 0.71 ±0.14 fb _1 . Since the statistical 
uncertainty is negligible, only the systematic one is quoted. 



7 D** production cross-sections 

The acceptance uncorrected distributions presented in Sec. [5] were given for all data taken in 2002-2006. 
The signals for D° and D° were summed up, and so were those for D* + and D* . In the following, the 
semi-inclusive differential cross-sections for production, determined for data from the year 2004 
only, will be obtained separately for D* + and D* . The acceptances, the integrated luminosity and the 
known branching ratio (2.6%) of D* to ^7171 are taken into account. At the end of this section, D* + 
and D* asymmetries will be shown for all 2002 to 2006 data, since integrated luminosity and also the 
acceptances cancel in these asymmetries to a good approximation. 

Figure [TT] displays the semi-inclusive differential cross-sections of D* + and D* events as a function of 
v,E,z and p\. The numerical values of the measured differential cross-sections are compiled in Table [I] 
These cross-sections are compared with the theoretical predictions obtained from the AROMA generator, 
which assumes cc production via photon gluon fusion and includes parton showers. The AROMA total 
cross-section is rescaled to the value of 1.9 nb measured by COMPASS, see below. They were calculated 
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Fig. 11: Semi-inclusive differential cross-sections for D* + and D* production as a function of (a) virtual 
photon energy v, (b) D° energy E, (c) fractional energy z and (d) squared transverse momentum p\. For 
all distributions, the red squares correspond to D* + and blue triangles to D* events (2004 data, D* sam- 
ple). The green circles are semi-inclusive differential cross-sections for D from the EMC experiment, 
see text. The curves represent AROMA predictions, dashed for D* and dotted for D* + . 



using the same program package parameters as those for the determination of acceptances. The p\ 
and z distributions are also compared with results published by the EMC Collaboration 20 years ago 
EH, based on 92 events, obtained with higher muon beam energy and a cut on Q 2 > 3 (GeV/c) 2 . EMC 
combined D° and D as within the statistical precision no differences were observed. In order to compare 
with the present data, their measured values and uncertainties are divided by a factor of 2. 

Good agreement is observed between the shapes of the measured distributions and the corresponding 
AROMA predictions. The distributions of D* + and D* as a function of v show that the points for D* 
are systematically higher than those for D* + . The effective threshold of D* + appears to be about 10 GeV 
higher than that of D* . The AROMA generator produces also somewhat more D*~ than D* + but the 
differences at threshold are far less pronounced. A similar feature can be observed for the z distribution. 
In the large-z region, that has a large contribution from low-v events, the cross-section of D* becomes 
significantly larger than that of D* + . The AROMA calculations predicts more D* than D* + as well, but 
the size of the effect is smaller. For the semi-inclusive differential cross-sections as a function of E and 
p\, no remarkable differences are observed between the shapes of the distributions of D* + and D* . 

The total cross-sections for D* + , D*~ and D*^ production are extracted by integration of the differential 
ones. The differences between the results from the integration over v, E, z and p\ (see Table [TJ are 
used to evaluate the systematic uncertainty of acceptance corrections. Using the RMS of the four results 
(from v, E, z and p\) one obtains a systematic contribution of 0.05 for both D* + and D* and 0.10 for 
the sum D**, i.e. at the level of the statistical uncertainty. In the ratio of D* + over D* the acceptances 
almost cancels. The values of the ratio vary between 0.77 and 0.81, with an average of 0.80 and a RMS 
of < 0.02, i.e. two to three times smaller than the statistical uncertainty of ~ 0.05. 
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Table 1: Semi-inclusive differential cross-sections for D* + and D* production as a function of (a) y* 
energy v, (b) D° energy E, (c) fractional energy z and (d) squared transverse momentum p\ of the D°. 
The central values and bin sizes of v and E are given in units of GeV, those of p\ in (GeV/c) 2 . The last 
two lines show the integrated cross-sections. Statistical uncertainties are given. 
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The final result for the D* meson production cross-section is then a(nN — > fl'D^X) = 1 .86 ±0.06 (stat) ± 
0.10 (sys) ±0.37 (luminosity) nb. The only cut applied is the energy window for the D° meson between 
20 GeV < E < 80 GeV in the laboratory frame, corresponding to 22 GeV < E < 86 GeV for the D* 
energy. 

For charm-anticharm production, AROMA gives a cross-section of 7.2 nb with 1.35 GeV chosen as 
the default charm quark mass. Using the common assumption of 0.6 D* mesons per charm event and 
accounting for the energy cut 20 GeV < E < 80 GeV, which reduces the number of charm Monte Carlo 
events by another factor of 0.6, the corresponding AROMA cross-section predicted for COMPASS is 2.6 
nb. Given the number of assumptions which underlie the AROMA default options (charm quark mass, 
fragmentation, no radiative corrections, leading order QCD apart from parton showers) the agreement 
with the above experimental result is considered to be good. 
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Fig. 12: Measured D* + and D* asymmetries for data (blue stars) and AROMA generator (red crosses) 
events as a function of X = v, E, z and p\ (coloured online). All 2002-2006 data are used. 



However, deviations from the AROMA predictions are observed in the data with respect to D* + and 
D* production. These may provide valuable insight into their production mechanisms. In a simple LO 
approach, assuming photon-gluon fusion with independent fragmentation of the charm and anti-charm 
quarks to be the relevant production mechanism, no differences should be observed between D* + and 
D* . Differences may occur for all processes where the quark content of the target nucleon matters. The 
quark content of D* mesons indicate that only D* may contain a valence quark from the target nucleon. 
Furthermore, instead of fragmenting into D* + the c quark, together with a diquark of the target nucleon, 
can hadronize into a charmed baryon, leading to associated production of e.g. D* A c . Thus the D* may 
result from a valence quark and/or associated production. If parton showers are included in AROMA the 
flavour dependent quark distribution functions of the nucleon come into play. Processes like associated 
production of D* A c lead to differences between kinematic distributions of D* + and D* . The same 
happens for processes where an initial quark in the nucleon absorbs the virtual photon and radiates a 
heavy gluon which then decays to cc, or where in the course of fragmentation the c quark picks up 
quarks from the nucleon. 

In order to provide statistically more precise information on the potentially interesting differences be- 
tween D* + and D* production, Fig. 12 shows particle-antiparticle asymmetries of the semi-inclusive 
cross-sections, 



Apr) 



do D ^{X)-dc D (X) 
dc D * + {X)+dc D> -{X) 



(6) 



as a function of X = v, E, z and p\ for both the D* sample and Monte Carlo events generated by 
AROMA. Here the full statistics of the years 2002-2006 is used. It is assumed that the acceptances for 
the two charge combinations are equal. In the previous section it was shown that for the year 2004 this 
is indeed approximately true. The numerical values of the measured asymmetries are given in Table |2| 
where only statistical uncertainties are shown, based on the assumption that acceptance cancels. A small 
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Table 2: Measured asymmetry A(X) as a function of X = v, E, z and p\. The central values and bin 
sizes of v and E in (a) and (b) are given in units of GeV, those of p\ in (d) in units of (GeV/c) 2 . 





(a) 




(b) 


v±Av/2 


A(v) 


E±AE/2 


A(E) 


28 + 4 


-0.130±0.415 


12.5 ±2.5 


-0.263 ±0.333 


36 + 4 


-0.610 + 0.098 


17.5±2.5 


-0.099 ±0.092 


44 + 4 


-0.272 ±0.082 


22.5 ±2.5 


-0.149 ±0.051 


52±4 


-0.207 ±0.059 


27.5 ±2.5 


-0.034 ±0.041 


60±4 


-0.032 ±0.051 


32.5 ±2.5 


-0.045 ±0.039 


68 ±4 


—0.174 ±0.047 


37.5 ±2.5 


—0.038 ±0.042 


76±4 


-0.123 ±0.044 


42.5 ±2.5 


-0.148 ±0.043 


84±4 


-0.078 ±0.046 


47.5 ±2.5 


-0.139 ±0.049 


92±4 


-0.1 16 ±0.053 


52.5 ±2.5 


-0.159 ±0.050 


100 ±4 


-0.109 ±0.048 


57.5±2.5 


-0.203 ±0.057 


108 ±4 


-0.027 ±0.057 


62.5 ±2.5 


-0.180±0.080 


U6±4 


-0.108 ±0.060 


67.5 ±2.5 


-0.103 ±0.109 


124 ±4 


-0.090 ±0.068 


72.5 ±2.5 


-0.235 ±0.185 


132±4 


-0.018 ±0.084 


77.5 ±2.5 


±0.074 ±0.231 


140 ±4 


±0.016±0.119 


82.5 ±2.5 


±0.593 ±0.604 


148 ±4 


-0.215 ±0.420 


87.5 ±2.5 


±0.288 ±2. 146 



(c) (d) 



z±Az/2 


Mz) 




p\ ± A/4 /2 


MPt) 


0.225 ±0.025 


-0.090 ±0.082 




0.10 + 0.10 


-0.094 ±0.046 


0.275 ±0.025 


+0.042 ±0.071 




0.31 ±0.11 


-0.105 ±0.046 


0.325 ±0.025 


+0.011 ±0.052 




0.56±0.13 


-0.151 ±0.045 


0.375 ±0.025 


-0.054 + 0.043 




0.84±0.15 


-0.082 ±0.045 


0.425 ±0.025 


+0.002 + 0.041 




1.16±0.17 


-0.083 ±0.045 


0.475 ±0.025 


-0.068 + 0.040 




1.53±0.20 


-0.076 ±0.044 


0.525 ±0.025 


-0.137 + 0.041 




1.96 ±0.23 


-0.148 ±0.049 


0.575 ±0.025 


-0.223 + 0.046 




2.45 ±0.26 


-0.192 ±0.058 


0.625 ± 0.025 


-0.193 + 0.052 




3.01 ±0.30 


-0.165±0.061 


0.675 ±0.025 


-0.260 + 0.059 




3.66±0.35 


-0.124 ±0.077 


0.725 ±0.025 


-0.479 + 0.078 




4.40 ±0.40 


-0.208±0.105 


0.775 ±0.025 


-0.609 + 0.099 




5.25 ±0.46 


-0.253±0.155 


0.825 ±0.025 


-0.939 + 0.213 




6.23 ±0.52 


+0.380 + 0.266 


0.875 ±0.025 


-0.926 + 0.419 




7.35 ±0.60 


-0.520 + 0.496 



cross-section assymmetry between D + and D production has been observed recenlty in a different 
energy range by the LHCb experiment |29l . 

As one can see from the figure, the measured asymmetry decreases significantly stronger than that pre- 
dicted by AROMA when v decreases below 40 GeV and/or when z increases above 0.6. The distributions 
shown as a function of v clearly exhibit different thresholds for D* + and D* production, which supports 
a stronger presence of mechanisms other than PGF with independent fragmentation. As a function of z, 
the most pronounced differences between D* + and D* are seen at large values of z, whereas at z values 
lower than 0.5 the production rates are nearly equal. Values of z larger than 0.5 indicate an asymmetric 
sharing of the energies between a D meson and its associated partner with opposite charm content. Since 
the cross-section of D* , which contains a down and an anti-charm quark, increases with increasing z 
stronger than that of D* + , this observation suggests processes where the anti-charm quark is fast and 
the charm-quark is slow. Here, a candidate process is again associated production of a D* along with 
a charmed baryon, i.e D* A c . Alternatively, since the D* may also contain a valence quark of the 
nucleon whereas the D* + does not, one may think of processes other than associated production, which 
involve valence quarks of the nucleon. 

Asymmetries between the production of D° and D° or D* + and D* were already observed in numer- 
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ous earlier experiments (see e.g. |22] for charm photoproduction and ll30"l - [36l for charm production by 
hadrons), although not as pronounced as in the present experiment that covers the region of virtual-photon 
energies from threshold up to 140 GeV. 

8 Summary and conclusions 

The observed total cross section of (1.9 ±0.4) nb for the production of D* + and D* mesons in in- 
elastic muon nucleon interactions at 160 GeV incident muon energy within the COMPASS acceptance 
(20 GeV < E < 80 GeV and 22 GeV < E < 86 GeV, for D° and D* respectively) lies within the range 
of values expected if the dominant process is photon-gluon fusion to open charm production. The total 
error is dominated by the uncertainty on the luminosity. 

The detailed comparison of the measured differential cross sections of D* production as a function of 
the variables v, E, z and p\ shows good agreement with those expected from the model underlying 
the AROMA generator used to produce the theoretical distributions. This is remarkable as most of the 
kinematic distributions of D mesons are quite different in shape compared to those of the background 
and the neighbouring X|(1430) resonance. 

The observed large asymmetries between D* + and D*~ production for v < 40 GeV and z > 0.6 can only 
partially be described by the model used in AROMA, which predicts differences of the same sign but 
of smaller magnitude. This indicates that the hadronization processes of charm and anti-charm quarks 
differ more significantly than expected or/and processes other than PGF contribute by a larger amount 
than assumed. 

The observed dependences of these differences on the kinematic variables, in particular on the photon 
energy v and the fractional energy z, suggest that associated production (e.g. DA C ) plays a dominant role 
at low photon energy. Also, D* production involving valence quarks of the nucleon may contribute to 
the observed asymmetries between D* + and D* production. 
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